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Singularities of Euler flow? Not out of the blue! 
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Does three-dimensional incompressible Euler flow with smooth initial conditions develop a sin- 
gularity with infinite vorticity after a finite time? This blowup problem is still open. After briefly 
reviewing what is known and pointing out some of the difficulties, we propose to tackle this is- 
sue for the class of flows having analytic initial data for which hypothetical real singularities are 
preceded by singularities at complex locations. We present some results concerning the nature 
of complex space singularities in two dimensions and propose a new strategy for the numeri- 
cal investigation of blowup. (A version of the paper with higher- quality flgures is available at 
http:/ /www. obs-nice.fr/etc7/complex.pdf) 



I. PHENOMENOLOGY OF BLOWUP 

According to Richardson's ideas on high Reynolds 
number three-dimensional turbulence, energy introduced 
at the scale i^, cascades down to a scale <C £o where it is 
dissipated. Consider the total time which is the sum of 
the eddy turnover times associated with all the interme- 
diate steps of the cascade. From standard phenomenol- 
ogy a la Kolmogorov 1941 (K41), the eddy turnover time 
varies as If we let the viscosity v, and thus 77, tend 
to zero, is the sum of an infinite convergent geometric 
series. Thus it takes a finite time for energy to cascade 
to infinitesimal scales, an observation first made by On- 
sager [1]. We also know that in the limit v ^ Q, the 
enstrophy, the mean square vorticity, goes to infinity as 
(to ensure a finite energy dissipation). 

From such observations, it is tempting to conjecture 
that ideal flow, the solution of the (incompressible) 3-D 
Euler equation 



dtv + V ■ Vv = — Vp, 
V-t> = 0, 



(1) 
(2) 



when initially regular, will spontaneously develop a sin- 
gularity in a finite time. 

This is of course incorrect: the kind of phenomenol- 
ogy assumed above is meant only to describe the (sta- 
tistically) steady state in which energy input and en- 
ergy dissipation balance each other; the inviscid — 0) 
initial-value problem is not within its scope. Another 
possible argument in favor of singularities has to do with 
the scaling properties of the high Reynolds number so- 
lutions (e.g. the spectrum). For the simpler case 
of the Burgers turbulence [2] , the scaling of spectra and 
structure functions is clearly rooted in the singularities 
(shocks) appearing in the solutions in the limit of van- 
ishing viscosity. It is however well-known that power-law 
behavior can be present without any singularities. An 



example is the Holtsmark process, that is any compo- 
nent of the electric or gravitational field produced at a 
given point by a set of charges or masses with an initial 
Poisson distribution in space and moving with uniform 
independent isotropic velocities (having, e.g., a Gaussian 
distribution). It is then easily shown (by adaptation of 
the technique used by Chandrasekhar [3]) that the cor- 
relation function is oc \t — 1'\~^. This power-law behavior 
comes from the algebraic distribution of the distances of 
closest approach to the point of measurement and not 
from singularities of individual realizations (which are 
actually analytic). 

There is yet another phenomenological argument, not 
requiring K41, which suggests finite-time blowup of the 
vorticity. Consider the equation for the vorticity u) = 
V A v for inviscid flow, written as 



DtU) = u) ■ Vt), 



(3) 



where Dt = dt + v - V denotes the Lagrangian derivative. 
Observe that Vv has the same dimensions as u} and can 
be related to it by an operator involving Poisson-type 
integrals. (For this use the fact that V^v = —V A w.) It 
is then tempting to predict that the solutions of (3) will 
behave as the solution of the scalar nonlinear equation 



(4) 



which blows up in a time l/s(0) when s(0) > 0. Actually, 
(4) is just the sort of equation one obtains in trying to find 
rigorous upper bounds to various norms when studying 
the well-posedness of the Euler problem. This is precisely 
why the well-posedness 'in the large' (i.e. for arbitrary 
i > 0) is an open problem in three dimensions. 

The evidence is that the solutions of the Euler equa- 
tion behave in a way much tamer than predicted by (4) 
because of a phenomenon known as "depletion of nonlin- 
earity" : when small-scale structures appear through the 
nonlinear dynamical evolution from smooth initial data. 
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they tend to display, at least locally, a much faster de- 
pendence on one particular spatial direction, so that the 
flow is to leading order one-dimensional [4]. An exam- 
ple in three dimensions are the vorticity pancakes seen in 
simulations, as illustrated in Figure 1. If the flow were 
exactly one-dimensional, the nonlinearity would vanish 
(as a consequence of the incompressibility condition). 




FIG. 1: Simulation of the formation of a vorticity pancake at 

four different times from Ref. [5] 

Depending on how strong this depletion is and also on 
how persistent it is, finite-time blowup may or may not 
occur [6, 7]. 

A review of the known mathematical results concern- 
ing the initial value problem for the Euler equation may 
be found in Ref. [8]. We just mention a few salient 
facts. In two dimensions, for sufficiently smooth initial 
data in a bounded domain (including periodic boundary 
conditions) or with sufficiently fast decay at large dis- 
tances. Holder continuity of the vorticity is preserved for 
all times. The key result was actually obtained in 1933 
[9, 10]. In three dimensions, for sufficiently smooth initial 
data, regularity is guaranteed only for a finite time. The 
first such result goes back to 1925 [11]. A very impor- 
tant result, established in the late eighties by Beale, Kato 
and Majda (BKAI) is that blowup, if it takes place, re- 
quires the time-integral of the supremum of the vorticity, 
and hence the vorticity itself, to become infinite [12]. As 
pointed out in Ref. [8] , the main stumbling block in trying 
to improve the existing 3-D regularity results beyond a 
finite time is our still very rudimentary understanding of 
the mathematics of depletion. It is worth stressing here 
that an even partial progress on depletion could play a 
crucial role in understanding regularity issues for the 3-D 
Navier-Stokes equation [13]. 

Numerical studies of 3-D blowup have been going on 
since at least the seventies with the development of spec- 
tral methods [14]. The advantage of spectral methods 
is that they allow the kind of very high accuracy which 
is desirable in investigating possible singular behavior. 
When spectral methods were able to achieve fairly high 
resolutions (256*^ or more grid points), evidence for de- 



pletion emerged, possibly of sufficient strength to prevent 
blowup [15] (see also Section III). It was however realized 
that singularities (or near singularities) of the 3-D Euler 
are highly localized in space and that the kind of uniform 
grid used in standard spectral methods is very waste- 
ful. A number of new methods were developed using 
nonuniform, e.g., adaptive meshing and often rather spe- 
cial initial conditions. Such studies have provided some 
evidence for blowup, based on extrapolating the behavior 
in time of the supremum of the vorticity. For the conve- 
nience of the reader we give a list of the key references 
(provided to us by R. Pelz and R. Kerr). It is not our 
purpose here to review such work. 

We shall argue in this paper that there exists a class 
of analytic flows for which numerical studies of blowup 
can probably be carried out with sufficient control to dis- 
tinguish genuine and spurious numerical blowup. In Sec- 
tion II, we recall some known facts about the Euler equa- 
tion in the complex domain. In Section III, we first recall 
how to detect precursors of blowup by tracing complex 
space singularities with a spectral method; then we pro- 
pose a new spectral adaptive approach capable, in prin- 
ciple, to resolve highly localized singularities. For this 
method it can be crucial to have information about the 
nature of complex space singularities. Sections IV and V 
report preliminary results concerning the singularities of 
the Euler equation in two space dimensions. Section VI 
presents concluding remarks. 



II. ANALYTICITY OF SOLUTIONS TO THE 
EULER EQUATION 

We are here interested in solution of the rf-dimensional 
Euler equation (2) with real analytic initial data vq, ex- 
tended into the complex domain, a question addressed 
also in Refs. [16, 17]. For simplicity, we assume "periodic 
boundary conditions" , that is space periodicity with pe- 
riod 1 in all d coordinates (although we use 27r when 
discussing numerical results) . The configuration space is 
thus = + iM."^, where T'^ is the d-dimensional peri- 
odicity torus. The analytic continuation of the solution 
to is denoted v{z,t), where z = x + iy (we do not 
complexify the time variable). 

We now give a brief survey of some key results for the 
Euler equation in the complex domain. Since this paper 
is intended mostly for a readership of fluid dynamicists, 
physicists and numerical analysts, we shall avoid using 
excessively formal mathematical language but, of course, 
distinguish clearly what is truly proven from what is just 
conjectured. 

The first results about analyticity of the solution to 
the Euler equation have been obtained, to the best of 

our knowledge, in the seventies. With periodic boundary 
conditions, analyticity assumed initially, is preserved for 
all time in two dimensions [18] and at least for a finite 
time in three dimensions [19, 20]. The simplest deriva- 
tions of such results are obtained by using Lagrangian 
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methods [21. 22]. One follows complex characteristics, 
the solutions of i = v{z, t), in order to control the width 
5{t) of the "analyticity strip" , that is the distance from 
the real domain of the nearest singularity in C'' for the 
solution at time t. 

What singles out two dimensions is that vorticity is 
conserved along the characteristics (as it is in the real 
domain). The key estimate used to prove all-time ana- 
lyticity is that when S{t) > s and < s < 1, one has 

\lmv{x + iy,t)\ < B,\\u}{-,t)\U\y\ In^, (5) 
||w(-,t)||s= sup \u>{x + iy,t)\, (6) 

\y\<s, xeT^ 

where Bg = Ci + C2 (ln(l/s))~"'^ and Ci and C2 arc pos- 
itive constants. This is a consequence of the Biot Savart 
law relating the velocity and the vorticity. Actually, (5) is 
mostly a reworking in the complex domain of an estimate 
obtained in Refs. [9, 10] for proving all-time regularity in 
the real domain. It follows from (5) and vorticity con- 
servation that, if initially (5(0) > s, at any later time 
the width of the analyticity strip has a double exponen- 
tial lower bound 5{f) > 5<=xp(B,,||t^o||.,t) Hence analyticity 
holds for all times, but complex singularities can get ar- 
bitrarily close to the real domain. 

In three dimensions, vorticity is not conserved since 
it can be stretched by velocity gradients. As a conse- 
quence, very poor control of d is available and its van- 
ishing after a finite time cannot be ruled out. Still, it 
is easy (but a bit technical) to show that, if at some 
time to the solution is analytic, it will stay so in a (pos- 
sibly small) time interval [to, to + Ti,[. Furthermore it 
was shown by Benachour and Bardos that, if 6{to) > s 
then, in this time interval, one has the square root lower 
bound S{t) > (s/2) (1 - (t - to)/^)^^^ [20-22]. One im- 
portant consequence, is that: In three dimensions with 
periodic boundary conditions and analytic initial condi- 
tions, analyticity is preserved as long as the velocity is 
continuously differentiable (C^) in the real domain [21]. 
The BKM theorem (cf. Section I) allows us to strength- 
ening this result: analyticity is actually preserved as long 
as the vorticity is finite. 

As long as 5 > the Euler equation can be written 
not only in the real domain but directly in the complex 
domain. In particular we shall find it useful to write it on 
"parareal domains" . By this we understand a domain of 
the form T''+iy for fixed y such that \y\ < 6, i.e. obtained 
from the real (periodic) domain T** by a fixed imaginary 
shift. In such a domain the velocity field is complex. 
By separating the velocity and the pressure into real and 
imaginary parts, v = v,: + ivi and p = p^ + ipi, we can 
rewrite the Euler equation on a parareal domain as 

dtV^ + ■ VVr — Vi ■ VVi = — Vpr. (7) 
dtVi + t)r • Vfi + Vi ■ VVr = -Vpi, (8) 

V • ■Ur = 0, V • -Ui = 0. (9) 



Note that eqs. (7-9) have some similarity to the magneto- 
hydrodynamics (MHD) equations, with and Vi playing 
the role of the velocity field and the magnetic field, re- 
spectively. The main differences are the following: (i) the 
presence of a pressure term in the second equation (there 
is no such term in the MHD induction equation); (ii) if 
we linearize the equation around a uniform Vi, we do not 
obtain Alfven-like waves but an instability whose growth 
rate is proportional to the wavenumber. The interpreta- 
tion of this instability is that the imaginary translation 
induced by an imaginary velocity corresponds to an expo- 
nential factor in Fourier space (rather than a phase factor 
which would be associated to a real translation) . We also 
observe that a real-imaginary decomposition similar to 
(7)-(9) has been used in Ref. [23] for the Navier-Stokes 
equation in the complex domain (in order to estimate the 
space analyticity radius of solutions in terms of and 
L°° norms of initial data). 

The proof of the local-in-time analyticity result of 
Ref. [19] is easily extended to the Euler equation in a 
parareal domain with analytic initial data in both two 
and three dimensions. All-time analyticity is now ruled 
out, even in two dimensions, since complex space singu- 
larities will generally cross after a finite time any parareal 
domain as they approach the real domain. 

A consequence of local-in-time analyticity is that the 
width of the analyticity strip d{t) cannot decrease dis- 
continuously in time. This is a consequence of the 
Benachour-Bardos square root lower bound on 5{t) given 
above. In particular if there is finite-time blowup, that is 
5{t) vanishes at some t^,, the vanishing takes place con- 
tinuously: real singularities do not come out of the blue 
(but out of the complex). This observation has led to 
the introduction of the method of tracing complex space 
singularities discussed in the next section. 

Note that similar results about finite-time analyticity 
and continuity of 5{t) can be obtained without periodic 
boundary conditions, provided the solutions decrease suf- 
ficiently fast at large distances [22]. However, if the so- 
lutions are not well behaved at infinity the results may 
be wrong. A counterexample, due to S. Childress and 
E. Spiegel (quoted in [24]), is given by 

, fx2+X3 X3+X1 Xi+X2\ 

v{x„X2, X3, t) = [-j^, -j^, -j-^ ) , (10) 

p{xi,X2,X3,t) = - ^^^^\\^2^ ■ (11) 

III. TRACING COMPLEX SINGULARITIES 

As suggested in Refs. [2, 25], when the initial data for 
the Euler equation are analytic (and periodic), it is pos- 
sible to trace the temporal behavior of the width of the 
analyticity strip 5{t) in order to obtain evidence for or 
against blowup. This takes advantage of the signature of 
complex space singularities in Fourier space. In one di- 
mension, it is well known that an analytic function having 
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isolated singularities in the complex plane has a Fourier 
transform whose modulus decreases at large wavenum- 
bers k as exp(— (5fc) (up to algebraic prefactors), where S 
is the distance from the real domain of the nearest com- 
plex space singularity [26, 27]. 

So far, nothing was known concerning the nature of 
complex space singularities of the multi-dimensional Eu- 
ler equation. In more than one dimension, singularities 
in the complex domain are never point-like. For example, 
it may be shown that for the d-dimensional (potential) 
inviscid Burgers equation with analytic initial data, the 
singularities in the complex domain (before shock for- 
mation) are on (d — l)-dimensional complex manifolds, 
a result which cannot be ruled out for the Euler equa- 
tion. In Fourier space, consider wavevectors of the form 
k = kk where k is unit vector of fixed direction. For 
fixed t, it may then be shown that, as fc — > oo, the mod- 
ulus of the Fourier transform decreases as exp(— (5(fe, t)k) 
where 6 now depends on the direction. The width of the 
analyticity strip at time t is then the minimum of 6{k, t) 
over all directions. 

The numerical tracing of complex singularities is easy 
if the Euler equation is integrated by a (pseudo-)spectral 
method [28] with enough spatial resolution to capture the 
exponential tails in the Fourier transforms. This method 
was applied for the first time to the three-dimensional 
Euler flow generated by the Taylor-Green initial condi- 
tions [15] 



vi= sm xi cos X2 cos , 
V2= — cosa;i sinx2 C0SX3, 
V3= 0. 



(12) 



The simplest is then to plot for various times the "energy 
spectrum" E{k,t), that is the angle-averaged squared 
modulus of the Fourier transform of the velocity. By 
a steepest descent argument one has (up to algebraic 
prefactors) E{k,t) oc exp(—S{t)k). This is illustrated in 
Figure 2, taken from Rcf. [15], giving for the first time 
evidence that the Taylor-Green vortex may actually not 
have any blowup since 5{t) appears to decrease exponen- 
tially in time. This behavior is observed reliably over 
an interval of time during which 6(t) decreases by about 
one decade. Later work, extending the simulations from 
a 256"^ to a 864^ grid, have confirmed this behavior over 
a range of one and a half decade in S (see Figure 3 taken 
from Ref. [5]). 

The Taylor-Green vortex may just happen not to be a 
good candidate for blowup. It is special in at least two 
ways. First, it has a lot of symmetry (used to simplify 
the computation). Second, its vortex lines have non- 
generic topology: they are closed instead of displaying 
the KAM disorder that is typical of the integral lines 
of three-dimensional divergenceless vector fields. Gen- 
eral periodic flow not displaying such pathologies has also 
been investigated in Ref. [5]. It seems to give exponen- 
tial shrinking of 6, but over just about half a decade in S. 
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FIG. 2: Spectral simulation of the inviseid Taylor Green vor- 
tex using 256'^ Fourier modes, (a) Evolution of the energy 
spectrum in lin log coordinates; from bottom to top: output 
from time t = 0.5 iu increments of 0.5. (b) Time-dependence 
of the width of the analyticity strip 5{t) in linear-log coordi- 
nates; the circles and plus signs correspond to 256^ and 128^ 
Fourier modes, respectively; the dashed line gives the thresh- 
old of reliability (from Ref. [15]). 




FIG. 3: Width of the analyticity strip 5{t) from a spectral 
simulation of a periodic flow with no particular symmetry. 
Dashed line data are for 200"^ modes and solid line data for 
256^ modes (from Ref. [5]). 



Two causes for this reduced range are clear: the limited 
resolution permitted ten years ago for flow without sym- 
metry (256'^) and a crossover phenomenon between two 
different small-scale structures localized at different spa- 
tial locations, resulting into two different regimes with a 
changeover around t = 1.05 (cf. Figure 3). 

If we could simulate general periodic flow in such a 
way as to observe the variation of S over several decades 
we would probably be able to get good evidence for or 
against blowup. One way to achieve this is just "pa- 
tience" . We can indeed expect to gain a little more than 
a factor four in spatial resolution every ten years: this 
requires an increase in CPU power of 4^ = 256, made 
possible by Moore's law. Another way is to renounce 
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spectral methods and switch to adaptive methods. Un- 
fortunately, such methods had so far only finite-order ac- 
curacy and could not be used to analytically continue 
the solutions into the complex domain, so that we must 
renounce measuring 5{t). 

We propose here a new strategy, the ■'spectral adap- 
tive" method which combines the highly localized refine- 
ment permitted by adaptive methods with full contact 
to the complex space structiue. This method, which is 
still in the testing phase, will only be briefly outlined 
here. The basic idea is to run a standard spectral simu- 
lation until the latest time ts when S{t) can be measured 
with very high accuracy and then to perform a "regulariz- 
ing analytic transformation" B on C'^ with the following 
properties: (i) it preserves T'^ globally and thus preserves 
periodicity, (ii) it maps the complex singularities of the 
solution at time ts away from (possibly to complex 
infinity). In the new coordinates resulting from the trans- 
formation B, the problem is still periodic and can again 
be integrated by a suitable spectral method (with new 
difficulties since the coefficients are now strongly non- 
uniform). The procedure can, in principle, be repeated 
several times. The spectral adaptive strategy is yet to 
be fully implemented. So far, we have performed tests 
in one space dimension on the Burgers equation. These 
have revealed that, in order to minimize errors, it is best 
to perform the regularizing transformation around the 
time when the round-off noise just disappears from the 
tail of the spatial Fourier transform of the solution, which 
is a function of the resolution and of the round-off level 
(cf. Figure 4). 




100 200 300 



wavenumber k ™^ 

FIG. 4: Evolution of the modulus of the Fourier transform 
for the solution of the one-dimensional Burgers equation with 
initial condition iio{x) = sin a;, for which the first real singu- 
larity is at t = 1. Around t = 0.82 round-off noise disappears 
from the resolved spectral range. A rather moderate reso- 
lution of 1024 modes is used in order to make the round-off 
noise fluctuations more visible. 

For one-dimensional equations, complex singularities 
are point-like and we know simple transformations which 



have the required properties. In higher dimensions, sin- 
gularities arc on extended objects and their nature is not 
well known. A possible candidate for the regularizing an- 
alytic transformation B could be the inverse Lagrangian 
map between time ts and the initial time (also called 
back-to-labels map [7]). This is definitely the case for 
the Burgers equation in any dimension, whose solution 
stays entire in Lagrangian coordinates if the initial con- 
dition is entire, that is analytic in the whole complex 
domain. However recent results for the two-dimensional 
incompressible Euler equation, the details of which will 
be published elsewhere, indicate that entire initial con- 
ditions develop complex-space singularities in both Eule- 
rian and Lagrangian coordinates. 

Note that the two-dimensional case is not just an aca- 
demic problem, as one might infer incorrectly from the 
proven fact that 2-D Euler fiow never blows up. One rea- 
son is that there is evidence that 2-D Euler flow is much 
tamer than predicted by the double exponential lower 
bound for S{t) given in Section 11. Indeed, spectral simu- 
lations indicate that 6{t) is actually decreasing exponen- 
tially [25] and that there is strongly depleted nonlinearity. 
Another reason is the existence of two-dimensional vari- 
ants of the Euler equation for which blowup has still not 
been ruled out, such as axisymmetric flow with swirl or 
2-D free convection [29-32]. 

To find suitable candidates for the analytic transforma- 
tion B in two dimensions we would like to know some- 
thing about the complex singularities of two-dimensional 
Euler flow with analytic initial conditions. This is the 
subject of the next section. 

IV. NUMERICAL STUDY OF COMPLEX 
SINGULARITIES FOR 2-D EULER FLOW 

Spectral simulations of 2D Euler flow with analytic ini- 
tial conditions give very strong evidence for the presence 
of complex space singularities through the presence of ex- 
ponential tails in the Fourier transforms of the solution. 
Can we find out something about the nature of such sin- 
gularities? One way is to analytically continue the so- 
lution at time t to complex locations, using the Fourier 
series. Suppose we have the Fourier representation 

t;(x,f)=^e*-%(f), (13) 

k 

we can just substitute z = x + iy for x and obtain v{z, t) 
as long as the series converges, that is for y < S{t). Al- 
ternatively, we can analytically continue the initial con- 
dition to the parareal domain -I- iy (a square in the 
2D case) and then integrate the parareal equations (7)- 
(9) from to t. Mathematically, the two procedures are 
equivalent since they differ only by an imaginary trans- 
lation, that is an overall exponential factor e~'°'^. Their 
numerical implementations in high-resolution spectral 
simulations may however not be equivalent, because of 
the presence of roundoff noise. At early times, when d{t) 
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can be quite large, the exponential falloff of the Fourier 
transform reaches roundoff noise level for wavenumbers 
which are much smaller than the maximum wavenum- 
ber permitted by the simulation (cf. Figure 5). If we 
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FIG. 5: Angle average, at a given time, of the modulus 
of the Fourier transform of velocity in lin-log coordinates. 
The exponential falloff reaches here the roundoff noise level 
for wavenumbers of the order of 350, while the maximum 
wavenumber femax determined by the 2/3 dealiasing rule is 
682 for a resolution of 2048^. 

multiply this by the exponential factor e"*"'^, roundoff 
noise may be tremendously amplified (for fc • y < 0). If 
we directly integrate in the parareal domain we can, at 
each time step, identify and set to zero the modes which 
are affected by roundoff, a procedure closely related to 
that used by Krasny for integration of the vortex sheet 
problem [33]. As t increases, these modes will shift to 
higher wavenumbers, eventually reaching the max;imum 
wavenumber permitted by the resolution. We are then 
effectively making better use of the available resolution. 
We call this procedure parareal integration with noise 
suppression. This method is rather crude and might be 
improved by taking into account the strong anisotropy 
present near a singularity. 

We have used this method to integrate the two- 
dimensional parareal Euler equation using from 512^ to 
2048^ modes and initial conditions which are trigonomet- 
ric polynomials in the space coordinates and thus entire 
functions. (As we shall see later, such initial conditions 
are also amenable to a short-time expansion.) The clean- 
est results are obtained when the initial condition vo{x) 
has only two modes 

vo{x) = {d2^po{x), -ditpoix)), (14) 
tpo{x) = cos(a;i) -|- cos(2a;2), (15) 

where ip{x,t) is the stream function and di and 82 are 
the derivatives with respect to a;i and X2- 

Figures 6 and 7 show the real part of the vorticity 
u! = — V^V-" at the latest time when the behavior in the 
parareal domain is still sufficiently smooth to be unaf- 
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FIG. 6: Snapshot of the contours of the real part of vorticity 
in the parareal domain y = (27r/64)(cos 7r/6, — sin 7r/6) just 
before the appearance of a singularity. 
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FIG. 7: Enlargement of Figure 6 showing the contour lines of 
the real part of vorticity in the neighborhood of the location 
where the singularity appears. 



fected by truncation, that is when the complex singular- 
ities are not too close to the parareal domain chosen. It 
is seen that the solution in the near-singular region is 
almost one dimensional and hence has strongly depleted 
nonlinearity. The structure of the region of high vorticity 
gradients suggests that the singular set is a smooth curve. 
A cross section of the vorticity in the near-singular region 
is shown in Figure 8; a blowup of the region of highest 
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vorticity (cf. Figure 9) shows about one decade of an 
approximately one-over-square root behavior (exponent 
— 1/2) of the vorticity when crossing the singular mani- 
fold. This implies some range of square-root behavior for 
the velocity. 

To check on errors due to truncation and filtering we 
changed the resolution from 512^ to 1024^ and 2048^ 
modes without changing any of the physical parameters 
(but, of course, adapting the time step to the spatial 
mesh) . The results diflFered by less than the filtering level, 
a good indication of reliability. We also tried to extend 
the scaling range at the higher resolutions by chosing a 
slightly later output time, so as to let the singularities 
move closer to the parareal plane. In fact, the scaling 
range did not increase appreciably, perhaps because of 
roundoff problems. We cannot therefore ascertain that 
the vorticity truly diverges with exponent —1/2 when 
approaching the singular set. 
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FIG. 9: Real part of the vorticity along the cut shown on 
Figure 8 in log-log coordinates. A range of algebraic behavior 
can be seen. 
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FIG. 8: Real and imaginary parts of vorticity along the cut 
represented as a dashed line on Figure 6. Both real and imag- 
inary parts of vorticity display quasi-singular behavior near 
s = —1.2 and s = 0. 



V. ASYMPTOTIC ANALYSIS OF 2-D 
COMPLEX SINGULARITIES 

In Ref. [34] the following result was established: the in- 
viscid one-dimensional Burgers equation with initial con- 
ditions which arc trigonometric polynomials in the space 
variable has, at short real times, square root branch point 
singularities for the velocity, which are within a distance 
d{t) oc \n{l/t) of the real domain. It was also briefly 
pointed out that the ln(l/t) law can probably be ex- 
tended to the 3-D Euler equation with initial conditions 
which are trigonometric polynomials in the space vari- 
ables (such as the Taylor Green flow). Further results 
on complex singularities for the Burgers equation may 
be found in Refs. [35, 36]. 

We now show that the same law applies to the 2-D Eu- 
ler equation at short times and, furthermore, we give a 



consistency argument shedding some light on the square 
root behavior reported in the previous section. For sim- 
plicity we shall limit ourselves to the two-mode initial 
condition (15) used in the previous section, although 
most of the arguments can be generalized to arbitrary 
trigonometric polynomials. 

We start with the 2-D Euler equation written in stream 
function formulation 

9tW = ./(V',VV), (16) 

where J(/, g) = d\fd2g — d^gdif- It is easily shown that 
(16) has a solution in the form of a temporal Taylor series 

^(a;,i) = ^^„(a;)r, (17) 

n>0 

where ipo is the initial condition and the ipnixY^ for n>l 
are defined recursively by 

VVn+l = ^ E •^(V'm,VVp)- (18) 

m+p=n 

For the two-mode initial condition tpo{x) = cos(a;i) -|- 
cos(2a;2), it is easily checked that tpn+i is a trigonomet- 
ric polynomial which, when written in terms of complex 
exponentials, involves only modes of the form Qvxi+2iqx2 
where p and q are signed integers such that |p| + |g| < 
n + 1. Each term can now be continued from real x to 
complex z = X + iy = (zi, Z2)- 

Trigonometric polynomials are instances of entire func- 
tions. Hence, initially, S = oo and, by continuity, 5{t) 
will be large at small real times. Singularities will thus 
be present only for suitably large |y|. For large posi- 
tive yi and 2/2, the dominant contributions to ipn{x + iy) 
have p < and 5 < 0. (For our special choice of initial 
conditions, symmetry arguments make it unnecessary to 
examine the other three sign quadrants in the y plane.) 
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When yi +00 and 2/2 — * +00, the dominant terms in 
tpn{z) are of the form 

n 

k=l 

where the complex coefficients satisfy suitable recur- 
sion relations, not needed here. If we formally keep only 
those dominant terms in (17), we obtain 

^{z,t) ~ ^ (e-«i +e-2*^=) 

n>Ofe=l 

Wc now let t ^ and simultaneously |y| oo in such a 
way that ie"*^^ and te"^*^^ stay finite and we find that (i) 
all the terms in the Taylor expansion stay finite and (ii) 
all the terms not included in (20) are subdominant. This 
observation leads us naturally to making the asymptotic 
ansatz ^p{z,t) ~ {l/t)F{z) where 



(5i, Z2) = {zi +i\nt, Z2 + {i/2)lnt). 



(21) 



Straightforward substitution into (16) leads to 

V2(-l + idi + {i/2%)F = J{F, V'F), (22) 

where the ovcrscript tilde means that the partial deriva- 
tives are taken with respect to the new variables. The 
initial condition (15) becomes now a boundary condition 



F{z) 



-2iZ2 



yi 



-00, 2/2 



-00. 



(23) 

If (22) with this boundary condition has a miique solu- 
tion possessing complex space singularities at a finite dis- 
tance from the real domain, then it follows immediately 
from the change of variables (21) that S{t) oc ln(l/i) as 
f — > in the original variables. We have recently checked 
numerically that, for the 2-D flow with initial condition 
given by (15), this scaling law holds over more than 8 
decades, up to f 0.1. Details will be pubhshed else- 
where together with numerical solutions of the asymp- 
totic equation (22). 

The numerical results of the previous section suggest 
that, for fixed t, the singularities in are located on one- 
dimensional complex manifolds, near which the velocity 
has some range of square root behavior. 

Using the short-time asymptotic equation (22), we 
show now that the observed behavior is consistent with 
the two-dimensional Euler equation. For this, assume 
that the (rescaled) stream function F{z) can be repre- 
sented, at least approximately, as 



F{z) = G{z) + $"(z) + h.o.t., 



(24) 



where G and $ are analytic functions, a is a non-integer 

exponent, $ (but not G) vanishes on the singular man- 
ifold iS and h.o.t. stands for "higher order terms", that 



is terms involving higher powers of The expansion 
(24) is somewhat reminiscent of the singular expansion 
used in Refs. [16, 17], except that (i) we do not a priori 
suppose that the singular manifold $ = moves with the 
flow and (ii) we shall determine the value of the exponent 
a. Substituting (24) into (22), we obtain 



V2(-1 + z^i + -d2)G - J{G, V^G) 



Ip; 



J{G,^)-i{di + -d2)^ 



a-3 



W(a - 1) J($, {di^y + (^2$)^)$^"-^ + h.o.t., (25) 

where 7 = a{a — l){a — 2). In (25) the most singu- 
lar term near S is that involving <I>"~^, which cannot be 
balanced by any other term. Hence, its coefficient must 
vanish, thereby constraining the functions G and $ to 
satisfy J(G, $) — i{di + ^82)^ = 0. For non-integer a, 
the term involving $^0-2 g^^g^ (jg^jjjjQ^ balanced by any 
other term unless it is actually analytic, that is 2q; — 2 is 
an integer. The smallest possible value is a = 3/2. The 
velocity has then exponent 1/2, that is a square root be- 
havior near S. This is here derived luider the assumption 
that $ vanishes linearly near the singular manifold. Oth- 
erwise a different scaling law may be obtained. 



VI. CONCLUSION 

It is clear that the issue of finite-time blowup is still 
open for initially smooth 3-D Euler flow (and also for 
3-D Navier Stokes flow). Wc have here proposed inves- 
tigating this issue within the more restricted class of ini- 
tially analytic flow, for which (hypothetical) real singu- 
larities are necessarily preceded by singularities in com- 
plex space. For this purpose wc propose a new spectral 
adaptive strategy which requires the tracking of complex 
singularities and the use, at suitable times, of a regulariz- 
ing map sending singularities too close to the real domain 
away from it. 

Finally, we should mention that the issue of experi- 
mental study of 3-D Euler blowup was discussed at the 
workshop. We noted the following: if a flow is started 
by standard methods such as a moving grid, the "initial 
conditions" will have considerable small-scale excitation 
(viscous boundary layers generated at solid boundaries) 
and would in fact become singular if this was extrapo- 
lated to zero viscosity. Alternative ways may be tried 
where the fiow is put in motion initially by body forces 
with no small-scale component, such as electromagnetic 
or (ultra) sonic stirring. 
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